Stability estimate for the Helmholtz equation with rapidly jumping coefficients

The goal of this paper is to investigate the stability of the Helmholtz equation in the high- frequency regime with non-smooth and rapidly oscillating coefficients on bounded domains. Existence and uniqueness of the problem can be proved using the unique continuation principle in Fredholm's alternative. However, this approach does not give directly a coefficient-explicit energy estimate. We present a new theoretical approach for the one-dimensional problem and find that for a new class of coefficients, including coefficients with an arbitrary number of discontinuities, the stability constant (i.e., the norm of the solution operator) is bounded by a term independent of the number of jumps. We emphasize that no periodicity of the coefficients is required. By selecting the wave speed function in a certain \resonant" way, we construct a class of oscillatory configurations, such that the stability constant grows exponentially in the frequency. This shows that our estimates are sharp.Comment: a) Added references, b) rewritten the introduction with a summary of the results/techniques of the paper, c) Corrected typo

Convergence analysis of energy conserving explicit local time-stepping methods for the wave equation

Local adaptivity and mesh refinement are key to the efficient simulation of wave phenomena in heterogeneous media or complex geometry. Locally refined meshes, however, dictate a small time-step everywhere with a crippling effect on any explicit time-marching method. In [18] a leap-frog (LF) based explicit local time-stepping (LTS) method was proposed, which overcomes the severe bottleneck due to a few small elements by taking small time-steps in the locally refined region and larger steps elsewhere. Here a rigorous convergence proof is presented for the fully-discrete LTS-LF method when combined with a standard conforming finite element method (FEM) in space. Numerical results further illustrate the usefulness of the LTS-LF Galerkin FEM in the presence of corner singularities

Wave number-Explicit Analysis for Galerkin Discretizations of Lossy Helmholtz Problems

We present a stability and convergence theory for the lossy Helmholtz equation and its Galerkin discretization. The boundary conditions are of Robin type. All estimates are explicit with respect to the real and imaginary part of the complex wave number $\zeta\in\mathbb{C}$, $\operatorname{Re}\zeta\geq0$, $\left\vert \zeta\right\vert \geq1$. For the extreme cases $\zeta \in\operatorname*{i}\mathbb{R}$ and $\zeta\in\mathbb{R}_{\geq0}$, the estimates coincide with the existing estimates in the literature and exhibit a seamless transition between these cases in the right complex half plane.Comment: 29 pages, 1 figur

A Posteriori Error Estimation for Highly Indefinite Helmholtz Problems

We develop a new analysis for residual-type aposteriori error estimation for a class of highly indefinite elliptic boundary value problems by considering the Helmholtz equation at high wavenumber as our model problem. We employ a classical conforming Galerkin discretization by using hp-finite elements. In [Convergence analysis for finite element discretizations of the Helmholtz equation with Dirichlet-to-Neumann boundary conditions, Math. Comp., 79 (2010), pp.1871-1914], Melenk and Sauter introduced an hp-finite element discretization which leads to a stable and pollution-free discretization of the Helmholtz equation under a mild resolution condition which requires only degrees of freedom, where denotes the spatial dimension. In the present paper, we will introduce an aposteriori error estimator for this problem and prove its reliability and efficiency. The constants in these estimates become independent of the, possibly, high wavenumber provided the aforementioned resolution condition for stability is satisfied. We emphasize that, by using the classical theory, the constants in the aposteriori estimates would be amplified by a factor

Numerical treatment of retarded boundary integral equations by sparse panel clustering

We consider the wave equation in a boundary integral formulation. The discretization in time is done by using convolution quadrature techniques and a Galerkin boundary element method for the spatial discretization. In a previous paper, we have introduced a sparse approximation of the system matrix by cut-off, in order to reduce the storage costs. In this paper, we extend this approach by introducing a panel clustering method to further reduce these cost

A posteriori error analysis for elliptic pdes on domains with complicated structures

Summary.: The discretisation of boundary value problems on complicated domains cannot resolve all geometric details such as small holes or pores. The model problem of this paper consists of a triangulated polygonal domain with holes of a size of the mesh-width at most and mixed boundary conditions for the Poisson equation. Reliable and efficient a posteriori error estimates are presented for a fully numerical discretisation with conforming piecewise affine finite elements. Emphasis is on technical difficulties with the numerical approximation of the domain and their influence on the constants in the reliability and efficiency estimate